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In order to understand the binding of four static quarks, flux distributions corresponding to these binding 
energies are studied in quenched SU(2) and also related to a model for the energies. The potential relevant to 
string breaking between two heavy-light mesons is measured in quenched SU(3) using stochastic estimates of light 
quark propagators with the Sheikholeslami-Wohlert action. 



Our goal is to understand multi-quark interac- 
tions from first principles. Since "multi-quark" 
means that the system can be decomposed into 
more than one color singlet, the simplest such sys- 
tem consists of four quarks. In this work we try to 
understand this simple case, relevant for meson- 
meson interactions, hoping that generalization to 
more complicated cases will be relatively straight- 
forward. 

We first give a brief introduction to a model 
that reproduces one hundred ground- and excited 
state energies of four static quarks [1] . These en- 
ergies have been calculated earlier in SU(2) gauge 
theory on rs 20 3 x 32 lattices with a « 0.1 fm - a 
continuum extrapolation has been performed [2] . 
The flux distribution producing four-quark bind- 
ing has been measured [3] and we attempt to re- 
late it to our model. 

Preliminary results are reported from a more 
realistic calculation concerning the potential be- 
tween two heavy-light mesons in quenched SU(3). 
The heavy quarks are taken to be static and 
the light quark propagators are calculated us- 
ing stochastic estimators with an 0(a) improved 
fermion action. This work also concerns string 
breaking. 

1. The model 

We like to think of the simplest multi-quark 
system, consisting of four quarks, as a bridge 
between a single hadron and interacting multi- 
hadron systems. Earlier, simulations of static 
four-quark systems have been performed for a 
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general set of geometries [2,4, references therein]. 
These quenched SU(2) calculations have pro- 
duced binding energies up to 120 MeV. In order 
to understand these results a phenomenological 
model based on two-body potentials and multi- 
quark interaction terms has been developed. 

The basis states of the model come from the 
three possible ways to pair four quarks into two 
mesons in SU(2). These pairing are taken to be 
in the ground state (A, B,C) or first excited state 
(A* ,B*,C*) of the two-body potential, giving six 
states (A, . . . , C*) in total. The basic equation of 
the model is 

[V - £(4)N] = 0, (1) 

where N contains overlaps of any two states 
(Nab — (A\B),...), V has the interactions 
(V AB = (A\V\B), . . .) and E(A) is the four-quark 
energy. The binding energy is defined as B(A) = 
E(4) — min(VAA, Vbb, Vcc), he. as the difference 
of the four-quark energy and the energy of the 
lowest-lying two-body pairing. 

A central element in the model is a multi- 
quark interaction term /, defined as {A\B) = 
(B\C) = (A\C) = -J/2. For small distances 
/ = 1 and for large / = 0. The potentials V 
have a perturbative one-gluon exchange prefac- 
tor and lattice two-body potentials vy , giving the 
form Vij = — | ^ Kj - • fjVij. This expression 
alone would lead to unphysical van der Waals 
forces, which are removed by the /. The in- 
teraction between basis states A and B is thus 
V AB = -{(V AA + V BB - V CC )- 

The factor / can be parameterized as / = 
exp(—b s kfS), where b s is the string tension and S 
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the minimum area bounded by the four quarks. A 
simple version of the model with just states A, B 
and / = 1 is reproduced by perturbation theory 
toO(g 4 )[5]. 

For quarks at the corners of a regular tetra- 
hedron binding increases with size, which sug- 
gests two-body excitations are relevant, as they 
are closer to the ground state for larger dis- 
tances. Therefore, we have included excited basis 
states in our model. This gives new /-factors 
(A*\B*) = ... = -r/2, (A*\B) = ... = -f a /2, 
the former of which is parametrized like / and the 
latter as f a — bsfgS exp(—b s k a S). The interac- 
tion matrix elements involving two-body excita- 
tions are parameterized in a manner similar to the 
ground state matrix elements, e.g. (A*\V\B*) = 
- S - [Vaa + V BB - V CC + c (VX A + V* B - VS C )}. 

When this model is fitted to 100 ground and ex- 
cited state four-quark energies it turns out that 
f c is always consistent with one, suggesting that 
two-body excitations interact in a perturbative 
manner. We are left with four independent pa- 
rameters k f = 1.5(1), k a = 0.55(3), /J! = 
0.51(2), Co = 3.9(3) that fit the energies with 
X 2 /d.o.f. w 1. The values of these parameters are 
stable when energies from smaller lattice spacings 
are used. If the excited basis states are left out 
we get x 2 ~ 3.2, which mostly comes from the 
maximally degenerate regular tetrahedra. 

2. Four-quark flux distribution 

The flux distribution is measured in quenched 
SU(2) using 

'(W(R,T)D^) - (W(R,T))(a^)' 



relate spatial sums over fields to energies via (3- 
functions such as b = dg/d\na = 2(S + U), e.g. 



JrK ' [ (W(R,T)) 

where W(R, T) is the Wilson loop, or more gener- 
ally, a combination of Wilson loops corresponding 
to the ground or excited state of the multi-quark 
system. Since extracting a signal is not easy we 
use link- integration on the Wilson loops. As in 
the calculation of energies, we utilize a variational 
approach in order to extract the ground or an ex- 
cited state of the system. The variational basis 
is formed of Wilson loops with various levels of 
fuzzing. 

We can check the accuracy of our flux mea- 
surement with lattice sum rules. These sum rules 
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^JS{£ x + £ y + £ z ) + U{B z + B y + B x 



Using our earlier two-quark results [6] for S, U 
this relation allows us to compare the sum over 
the distribution to a much more accurately mea- 
sured energy. This is especially convenient when 
the binding is considered; E is a lattice artefact 
from quark self-energies, and is removed for the 
binding energy. 

Three-dimensional 
pictures of our results (Fig. 1) can be viewed at 
www . physics . helsinki . f i/~ppennane/pics/. 
The binding distribution is found to have approx- 
imately constant height in between the quarks, 
whereas the first excited state has a "cloverleaf" 
symmetry with sign changes. Replacing energies 
with flux distributions in our model without exci- 
tated basis states gives kf < 1, which agrees with 
fits to energies [3]. When excited basis states are 
included the ground-state interaction has shorter 
range; kf > 1, suggesting that ground state po- 
tentials are important at small distances and ex- 
cited ones at larger distances. 
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Figure 1. The distribution of energy around four 
static quarks in a square with side R — 4 corre- 
sponding to a) total energy b) binding energy c) 
first excited state binding energy. 
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3. Two heavy-light mesons 

We calculate the potential between two heavy- 
light mesons as a function of the heavy quark sep- 
aration R in quenched SU(3) at j3 = 5.7. The 
heavy quarks are taken to be static and the light 
ones have approximately the s mass (Mps/Mv = 
0.65) with the propagators obtained as stochas- 
tic pseudofcrmionic estimates with maximal vari- 
ance reduction [7] . The principal advantage of the 
stochastic estimates is the ability to get propaga- 
tors to and from most points on the lattice; hence 
the nickname "all-to-all" . 

At R = the meson-antimeson (BB) system 
should have a mass similar to a pion, while the 
BB system should resemble a heavy-light baryon; 
the former is indeed observed. The preliminary 
results in Fig. 2 do not show the one-7r exchange 
effect of dcuson models which leads to differences 
in the BB and BB potentials at large distances. 
As expected, at small distances the BB potential 
is more attractive. 

As discussed in other papers in these proceed- 
ings [8], the heavy-light meson-antimeson sys- 
tem is relevant to string breaking. Unfortunately 
calculations of this system at large enough dis- 
tances are very hard, but the all-to-all propagator 
scheme offers a way to overcome the difficulties. 
In our case the BB and static QQ energies be- 
come equal at R 1.3 fm [9]; this is reached 
by our ongoing calculations on a 16 3 x 24 lat- 
tice. A variational operator mixing analysis [10] 
similar to that presented by Philipscn in a Higgs 
model is then performed. By running our code 
on UKQCD unquenched configurations we hope 
to shed light on the differences due to sea quarks 
in the mixing region. 
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Figure 2. Preliminary results for potential be- 
tween two B mesons or a B and B - the data 
marked with a "2" includes light-quark exchange. 
Here aM B = 0.863(1). 
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